A Froggatt-Nielsen Model for 
Leptophilic Scalar Dark Matter Decay 

Christopher D. Caroneu and Reinard Primulandql 

High Energy Theory Group, Department of Physics, 
College of William and Mary, Williamsburg, VA 23187-8795 

(Dated: May 2011) 

Abstract 

We construct a model of decaying, TeV-scale scalar dark matter motivated by data from the 
PAMELA and Fermi-LAT experiments. By introducing an appropriate Abelian discrete symmetry 
and an intermediate scale of vector-like states that are responsible for generating lepton Yukawa 
couplings, we show that Planck-suppressed corrections may lead to decaying dark matter that is 
leptophilic and has the desired lifetime. The dark matter candidate decays primarily to lepton/anti- 
lepton pairs, and at a subleading rate to final states with a lepton, anti-lepton and standard model 
Higgs boson. We show that the model can reproduce the observed positron flux and positron 
fraction while remaining consistent with the bounds on the cosmic ray antiproton flux. 
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I. INTRODUCTION 



A number of earth-, balloon-, and satellite-based experiments have observed anomalies in 
the spectra of cosmic ray electrons and positrons. Fermi-LAT l| and H.E.S.S. ^ have mea- 
sured an excess in the flux of electrons and positrons up to, and beyond 1 TeV, respectively. 
PAMELA [sl, which is sensitive to electrons and positrons up to a few hundred GeV in 
energy, detects an upturn in the positron fraction beginning around 7 GeV, in disagreement 
with the expected decline from secondary production mechanisms. Recent measurements 
at Fermi-LAT support this result kl. In contrast, current experiments observe no excess in 

n n 

the proton or antiproton flux [5|. Although astrophysical explanations are possible [6[, these 
observations can be explained if the data includes a contribution from the decays of unsta- 
ble dark matter particles that populate the galactic halo 0]. The dark matter candidate 
must be TeV-scale in mass, have a lifetime of order 10^^ seconds, and decay preferentially 
to leptons. A number of scenarios have been proposed to explain the desired dark matter 
lifetime and decay properties j8-12|. 

To be more quantitative, consider a scalar dark matter candidate x which (after the 
breaking of all relevant gauge symmetries) has an effective coupling g^ff to some standard 
model fermion / given by QeffXfLfR + h.c. To obtain a lifetime of 10^^ seconds, one finds 
Qeff ~ 10~^^ if ~ 3 TeV. From the perspective of naturalness, the origin of such a 
small dimensionless number requires an explanation. One possibility is that physics near 
the dark matter mass scale is entirely responsible for the appearance of a small number, 
as is the case in models where a global symmetry, that would otherwise stabilize the dark 
matter candidate, is broken by instanton effects of a new non-Abelian gauge group Go- A 
leptophilic model of fermionic dark matter along these lines was presented in Ref. js]: the new 
gauge group is broken not far above the dark matter mass scale and the effective coupling 
is exponentially suppressed, geff oc exp(— 167r^/(7|,), where qd is the Gd gauge coupling. 



(An example o: 
found in Ref. 



" a supersymmetric model with anomaly-induced dark matter decays can be 
9[.) On the other hand, the appearance of a small effective coupling can 
arise if the breaking of the stabilizing symmetry is communicated to the dark matter via 
higher-dimension operators suppressed by some high scale M. Then it is possible that Qe/f is 
suppressed by {m^/My, for some power p; it is well known that for ~ (9(1) TeV and p = 
2, the correct lifetime can be obtained for M ~ (9(10^^) GeV, remarkably coincident with the 
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grand unification (GUT) scale in models with TeV-scale supersymmetry (SUSY) 10|]. If the 
LHC fails to find SUSY in the coming years, however, then the association of 10^^ GeV with 
a fundamental mass scale will no longer be strongly preferred. Exploring other alternatives 
is well motivated from this perspective and, in any event, may provide valuable insight into 
the range of possible decaying dark matter scenarios. 

The very naive estimate for g^ff discussed above presumes that the result is determined 
by a TeV-scale dark matter mass m^, a single high scale M and no small dimensionless 
factors. Given these assumption, the choice M = M*, where M* = 2 x 10^^ GeV is the 
reduced Planck mass, would not be viable: the dark matter decay rate is much too large for 
p = 1 {i.e., there would be no dark matter left at the present epoch) and is much too small for 
p = 2 [i.e., there would not be enough events to explain the cosmic ray excess). However, 
Planck-suppressed effects arise so generically that we should be careful not to discount them 
too quickly. What we show in the present paper is that Planck-suppressed operators can 
lead to the desired dark matter lifetime if they correct new physics at an intermediate scale. 
In the model that we present, this is the scale at which Yukawa couplings of the standard 
model charged leptons are generated via the integrating out of vector-like states. This sector 



will have the structure of a Froggatt-Nielsen model [13|] : an Abelian discrete symmetry will 
restrict the couplings of the standard model leptons and the vector-like states, but will be 
spontaneously broken by the vacuum expectation values (vevs) of a set of scalar fields {</>}. 
Integrating out the heavy states will not only lead to the standard model charged lepton 
Yukawa couplings, but also to dark matter couplings that are naturally leptophilic and lead 
to dark matter decay. Aside from setting the overall scale of the charged lepton masses, the 
symmetry structure of our model will not restrict the detailed textures of the standard model 
Yukawa matrices. This feature is not automatic; symmetries introduced to guarantee dark 
matter leptophilia may also make it difficult to obtain the correct lepton mass matrices, at 
least without additional theoretical assumptions (for example, the addition of electroweak 
Higgs triplets, as in the model of Ref. 12] )• Our framework is free of such complications 
and is compatible, in principle, with many possible extensions that might address the full 
flavor structure of the standard model. 

Our paper is organized as follows. In the next section, we present a model that illustrates 
our proposal. In Section 3, we compute the predicted flux, $(e^), and the positron 
fraction $(e"'")/[$(e"'') + $(e~)] for some points in the parameter space of our model and 
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compare our results to the relevant cosmic ray data. It is worth noting that this analysis has 
applicability to any model that leads to similar dark matter decay operators. In Section 4, 
we comment on the relic density and dark matter direct detection in our example model. In 
Section 5, we summarize our conclusions. 



II. A MODEL 



We assume that the right-handed charged leptons of the standard model, cr, and four 
sets of heavy vector-like charged leptons are constrained by the discrete symmetry 

G = Zp X , (2.1) 

with p and q to be determined shortly. We assume that the vector-like leptons have the 
same electroweak quantum numbers as cr 

E^S^Ef^CR, (2 = 1... 4). (2.2) 

All the fields shown are assumed to be triplets in generation space, with their generation in- 
dices suppressed. Under the discrete symmetry, the fields in Eq. (12. 2p are taken to transform 

as 

eR->w"'^eR, (2.3) 

(2 = 1... 4). (2.4) 

We will take u and rj to be elements of Zp and Z^, respectively, with = 1 and rj'^ = 1. In 
addition, we assume the presence of a heavy right-handed neutrino, ur, that is a singlet under 
G. We note that the fields that are charged under G do not transform under any of the non- 
Abelian standard model gauge group factors, so that G satisfies the consistency conditions 
of a discrete gauge symmetry in the low-energy theory |l^; such discrete symmetries are 
not violated by quantum gravitational effects^. The Yukawa couplings of the standard 



^ The consistency conditions require that anomahes involving the non-Abelian gauge groups that are hnear 



in a continuous group that embeds G must vanish, as is automatic above. Ref. |14| indicates that no 
rigorous proof exists that the cancellation of the linear gravitational anomalies is a necessary condition 
for the consistency of the low-energy theory. Nonetheless, such a cancellation can be achieved here by 
including a singlet, left-handed fermion, A'^^. that transforms in the same way as under G. For the 
choice p = 8, adopted later in this section, Nl can develop a Majorana mass somewhat below M, and 
decay rapidly to lighter states via Planck-suppressed operators. Including such a state does not affect the 
phenomenology of the model otherwise. 
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model charged leptons arise when the symmetry G is spontaneously broken and the vector- 
like leptons are integrated out of the theory. Symmetry breaking is accomplished via the 
vacuum expectation values of two scalar fields 0e and (pOi which transform as 

<t>D Tl<t>D- (2.5) 

The following renormalizable Lagrangian terms involving the charged lepton fields are al- 
lowed by the discrete symmetry: 

i=l 

4 

+ Y^M^'^^f^R + h-c. (2.6) 

1=1 

While it is not our goal to produce a theory of flavor, we note that the terms in Eq. (12.61) are 
of the type one expects in flavor models based on the Froggatt-Nielsen mechanism. Hence, 
integrating out the E flelds leads to a higher-dimension operator 

CD^LLH<p%eR+ h.c, (2.7) 

which provides an origin for the charged lepton Yukawa couplings. Choosing {(j)E)/M ~ 0.3 
gives the correct scale for the tau lepton Yukawa coupling; the smaller, electron and muon 
Yukawa couplings may be accommodated by suitable choices of the undetermined couplings 
in Eq. (12. 6p . One might imagine that the remaining Yukawa hierarchies could be arranged 
by the imposition of additional symmetries, though we will not explore that possibility here. 
We now introduce our dark matter candidate x? a complex scalar fleld that transforms 

as 

X-^w^X and x^V~'^X (2-8) 

under xZq. We assume that all the nonvanishing powers of u and rj shown in Eqs. (12.31) . 
(12. 4p and (12. 8p are nontrivial, which requires that p > 4 and q > 2. Then, there are no renor- 
malizable interactions involving a single x fi^ld (or its conjugate) and two fermionic flelds 
that could lead to dark matter decay. However, non-renormalizable, Planck-suppressed oper- 
ator provide the desired effect. The lowest-order, Planck-suppressed correction to Eqs. (12. 6p 
that involves a single x fi^ld is the unique dimension-six operator 

AC^ = ±^xE?4>leR+ h.c. (2.9) 



Including Eq. (12. 9p and again integrating out the heavy, vector-like states, one obtains a 
new higher- dimension operator, 

(Idecay = XLlHcr + h.C, (2.10) 

which leads to dark matter decay. For ~ 3 TeV (compatible qualitatively with fits to 
the PAMELA and Fermi-LAT data), a hfetime of 10^^ seconds is obtained when 

~ 1 X 10-^^ (2.11) 

For our operator expansion to be sensible, we require {(po) < M; however, we also do not 
want a prohferation of wildly dissimilar physical scales, if this can be avoided. Interestingly, 
if we choose M to be the geometric mean of (H) and M*, one finds 

M = 2 X 10^° GeV, {(^e) = 0.3 M, (0z)) = 0.1 M , (2.12) 

which meets our aesthetic requirements. Standard model quark and neutral lepton masses 
are unaffected by the discrete symmetry of our model, by construction. Light neutrino 
masses arise via a conventional see-saw mechanism, and it is possible to obtain a right-handed 
neutrino mass scale Mr ^ M, so that all the heavy leptons appear at a comparable scale. 
Assuming that the largest neutrino squared mass is comparable to Amgg = 2.43 x 10~^ eV^, 



as suggested by atmospheric neutrino oscillations [15|, then this possibility is obtained if 
the overall scale of the Yukawa coupling matrix that appears in the neutrino Dirac mass 
term is of the same order as the charm quark Yukawa coupling. This scenario is depicted 
in Fig. [H In this case, the theory is characterized by three fundamental scales: the Planck 
scale, an intermediate scale (associated with charged lepton flavor and right-handed neutrino 
masses), and the TeV-scale. Symmetry-breaking vevs appear within a factor of < 10 below 
the latter two. Of course, the right-handed neutrino scale need not be linked with the scale 
at which the charged lepton Yukawa couplings are generated; this is simply one of many 
viable possibilities that depend on choices of the free parameters of the model. 

Finally, we return to the discrete symmetry group G = Zp x Z^. We have noted that the 
structure of the theory that we have described is obtained for p > A and q > 2, but this 
does not take into account an important additional constraint: there must be no Planck- 
suppressed operators involving couplings between the various scalar fields in the theory that 
can lead to other dark matter decay channels that are either (i) too fast or (ii) too hadronic. 
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FIG. 1: A possible choice for the mass scales in the theory. Symmetry breaking vevs appear within 
approximately an order of magnitude of the lower two scales. 

For example, the choice p = 5 and g = 3, allows the renormalizable G-invariant operator 
which leads to mixing, for example, between the x ^-^d 0s fields; the latter couples 
to two standard model leptons via the operator in Eq. (12. 7p . leading to a disastrously large 
decay rate. We find that all unwanted operators are sufficiently suppressed if we take p = 8 
and g = 4, that is 

Gi = ZsxZ^. (2.13) 

The lowest-order combination of scalar fields that is invariant under Gj, as well as the 
standard model gauge group, is 

(2.14) 



M3 



Suppression by three factors of the Planck scale is more than sufficient to suppress any 
operators that are generated when the (pE and </)£, fields are integrated out of the theory, or 
that may be constructed from products of Eq. (I2.14p with any Gj-singlet, gauge-invariant 
combination of standard model fields. It is straightforward to confirm that the alternative 
choice 

Gn = Z8XZ,, (2.15) 

is also viable, by similar arguments. The difference between the symmetry groups Gj and 
Gjj is that the former allows two types of dark matter mass terms: + h.c. and x^X- 
This leads to a mass splitting between the two real scalar components of so that the 
lighter is the dark matter candidate. The choice Gn forbids the mass terms, so that 
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the dark matter consists of particles and anti-particles associated with the original complex 
scalar field. We note that in this theory, the renormalizable interactions involving x have an 
accidental U(l);^ global symmetry which would lead to dark matter stability in the absence 
of the Planck-suppressed effects. The analysis that we present in the following sections is 
somewhat simplified by the choice of G//, which we adopt henceforth. 



III. COSMIC RAY SPECTRA 



In this section, we investigate the cosmic ray and proton/antiproton spectra of our 



model. Our treatment of cosmic ray propagation follows that of Ref. [16|. We show that 
model parameters may be chosen to accommodate the positron excess and the rising electron- 
positron fiux observed by the PAMELA and Fermi-LAT experiments, respectively. 

In Eq. fl2.10p . we identified the operator responsible for dark matter decays. More explic- 
itly, this operator may be written 

Cdecay = (^ij^^2 X^lHc^r + h.C, (3.1) 

where i and j are generation indices, and Cjj represents unknown order-one coefficients. 
Different choices for the couplings Cij will lead, in principle, to different cosmic ray spectra. 
To simplify the analysis, we focus on two possibilities: In the lepton mass eigenstate basis, 
the fermions appearing in the decay operators are either [i) muons exclusively, or (ii) taus 
exclusively. We will find that either of these choices is consistent with the data, even though 
we have not fully exploited the parametric freedom available in the Cjj. This is sufficient to 
demonstrate the viability of our model. The remaining factors in the operator coefficient are 
chosen to obtain the desired dark matter lifetime, as we discussed in the previous section. 
In unitary gauge, the operator (13. ip can be be expanded 

Cdecay = -^gij{Venj + h) X^^L (^R + ^.C, (3.2) 

where h is the standard model Higgs field, which we will assume has a mass of 117 GeV, 
Vew = 246 GeV, and gij = Cij{(l)D)'^ / {MM^). The term proportional to the Higgs vev leads 
to the two-body decay x ~^ ^ for i = fi or t, while the remaining term contributes to 
X — > i'^i^h. We take both of these decay channels into account in our numerical analysis. 
The final state particles in these primary decays will subsequently decay. The electrons. 
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positrons, protons and antiprotons that are produced must be added to expected astrophys- 
ical backgrounds to predict the spectra at experiments hke PAMELA and Fermi-LAT. 

Electrons and positrons that are produced in dark matter decays must propagate through 
the Milky Way before reaching the Earth. In order to determine the observed fluxes, one 
must model this propagation. The transport equation for electron and positrons is given by 

= V ■ [K{E,f)VU] + ^ [b{E,f)f,±] + Qe±(^,r), (3.3) 

where ff>±{E,'r,t) is the number density of electron or positrons per unit energy, K{E,f) is 
the diffusion coefficient and b{E, r) is the energy loss rate. We assume the MED propagation 



model described in Ref. 



171 ]. The diffusion coefficient and the energy loss rate are assumed 



to be spatially constant throughout the diffusion zone and are given by 

K{E, r) = O.OlUe^-™ kpcVMyr (3.4) 

and 

6(E, f) = lO^^V GeV/s , (3.5) 
where e = E/1 GeV. The last term in Eq. fl3.3p is the source term given by 

where is the dark matter mass and is the dark matter lifetime. In models like ours, 
where the dark matter can decay via more than one channel, the energy spectrum dN/dE 
is given by 

dN_^V,(dN\ 

Je-1^t\je)^^ (3.7) 
where Fj/F is the branching fraction and (dN/dE)i is the electron-positron energy spectrum 



of the i decay channel. We use PYTHIA [18| to determine the {dN/dE)i. For the dark 



matter density, p(r), we adopt the spherically symmetric Navarro- Frenk- White halo density 
profile Q 

(r/r,)[l + (r/r,)]2 

with po — 0-26 GeV/cm^ and ~ 20 kpc. The solutions to the transport equation are 
subject to the boundary condition /e± = at the edge of the diffusion zone, a cylinder of 
half-height L = 4 kpc and radius -R = 20 kpc measured from the galactic center. 
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The solution of the transport equation can be written 



ME}-J^jyE'GME,E')'-^, (3.9) 



where Ge±{E,E') is a Green's function, whose exphcit form can be found in Ref. 20|]. The 
interstellar flux then follows immediately from 

^fi' = ^fe4E). (3.10) 



16|: 



We adopt a parameterization of the interstellar background fluxes given in Ref. 

<^iE) = ( ^^°L,^, J GeV-m-s-'sr-', (3.11) 



1 + 0.224e2-93 



$^^/(^) = ( ^^-^^ '''' + 2A.0e~' A GeV-'m~\-hT-\ (3.12) 
^+ ^ ^ Vl + 0.0002e5-63 J ^ ' 

Finally, the flux at the top of the earth's atmosphere, ^^±^, is corrected by solar modulation 



effects 



16| 



<^IF\EroA) = ^<l>iUE,s) , (3.13) 

-C/r 



El 

where Eis = Etoa + |e|0, and |e|0 = 550 MeV. Eis and Etoa are the energy of 
positron/electron at the heliospheric boundary and at the top of atmosphere, respectively. 
The total electron and positron flux is determined by 

$*°*(E) = ^^^^\E) + $f/^(^) + k^^^.\E) + $?(^), (3.14) 

where is a free parameter that determines the normalization of the primary electron flux 
background. The positron excess is given by 

$«/^(E) + <l>'^^(E) , , 

P^E) = \L^Ef' • (3-15) 

The results of our analysis are presented in Figs. [2] and [3l In the case where the dark 
matter decays only to fi~^fi~ and fi~^fi~h, we find good agreement with the data for = 
1.8 X 10^^ s and = 2.5 TeV. In this case, the branching fraction to the two-body decay 
mode is 90.2%. In the case where the decay is to r+r^ and r^r^h only, our best results 
are obtained for = 9.0 x 10^^ s and = 5 TeV, corresponding to a two-body branching 
fraction of 69.6%. In all these results, the background electron flux parameter k is set to 



0.88, following Ref. |20|. 
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FIG. 2: Left panel: The positron excess for dark matter decaying into fJ.'^fJ," and The dark 

matter mass is 2.5 TeV and Hfetime 1.8 x 10^^ s; the branching fraction to the two-body decay mode 
is 90.2%. The dashed Hne represents the background and the sohd hne represents the background 
plus dark matter signal. Data from the following experiments are shown: PAMELA Q] (solid 



dots), HEAT 



2l| 



AMS-01 



M (V), and CAPRICE 



23( 1 (A). Right panel: The corresponding 
graph for the total electron and positron flux. Data from the following exp eriments are shown: 
Fermi-LAT ^ (solid dots), HESS ^ (v), PPB-BETS [20] (o), HEAT [23] (A). 



Since the dark matter decays in our model include the production of standard model 
Higgs bosons in the final state, it is worthwhile to check that subsequent Higgs decays do 
not lead to an excess of cosmic ray antiprotons, in conflict with the experimental data. This 
will not be the case at our two benchmark parameter choices since the branching fraction to 
the three-body decay mode is suppressed compared to the two-body mode. The procedure 
for computing the cosmic ray antiproton flux is similar to that of the cosmic ray electrons 
and positrons. The transport equation for antiproton propagation within the Milky Way is 
given by 

= V ■ [k{T, r)Vfp - K.(r1/p-] + Qp-(T, r) (3.16) 

where T is the antiproton kinetic energy, Vc{r) is the convection velocity, and the source 
term Qp has the same form as Eq. (13. 6p . As in the case of propagation, the antiproton 
number density can be expressed in terms of a Green's function 

1 



fp{T) 



MxTx JO 



Tmax flN-(T') 



(3.17) 



where Gp{T,T') can be found in Ref. 20|. The relation between the antiproton number 
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FIG. 3: Left panel: The positron excess for dark matter decaying into T~r"'" and T~T~h. The dark 
matter mass is 5.0 TeV and Hfetime 9.0 x 10^^ s; the branching fraction to the two-body decay 
mode is 69.6% . Right panel: The corresponding graph for the total electron and positron flux. 

density and the interstellar flux of antiproton is given by 



$f -^(T) = |-/,(T) , (3.18) 



where v is the antiproton velocity. We also take account the solar modulation effect on the 
antiproton flux at the top of atmosphere, which is given by 

$r^(T™^) = ( yT^^f^^ l H'iTisl (3.19) 



where Tjs and Ttoa are the antiproton kinetic energies at the heliospheric boundary and at 
the top of atmosphere, respectively, with Tjs = Ttoa + For the proton and antiproton 



flux, we adopt the background given in Ref. 



28| 



Again assuming the MED propagation model Ref. [17j, we compute the antiproton flux 
and the antiproton to proton ratio for dark matter decays to fi~fi^ and fi"fi~^h, shown in 
Fig. m and for decays to r^r^ and r^r^h, shown in Fig. [5l We see that in both cases, the 
antiproton excess above the predicted background curves is small and consistent with the 
data shown from a variety of experiments. 



IV. RELIC DENSITY AND DIRECT DETECTION 



In this section, we show that the model we have presented can provide the correct dark 
matter relic density while remaining consistent with the direct detections bounds. The part 

12 
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FIG. 4: Left panel: The antiproton flux for dark matter decaying into /i^/i^ and fi^fM^h. The dark 
matter mass is 2.5 TeV and Hfetime 1.8 x 10^^ s; the branching fraction to the two-body decay mode 
is 90.2%. The dashed hne represents the background and the soHd hne represents the background 
plus dark matter signal. Data from the following experiments are shown: PAMELA 

Isi] 



dots), WiZard/CAPRICE 



2g| (solid 



301 1 (o), and BESS [31] (A). Right panel: The corresponding graph for 



the antiproton to proton ratio. Data from the following experiments are shown: PAMELA 



(solid dots), IMAX 



3| (*), CAPRICE 



30( (o) and BESS 



29| 



3l| (A). 




FIG. 5: Left panel: The antiproton flux for dark matter decaying into t^t^ and T~T~h. The dark 
matter mass is 5.0 TeV and lifetime 9.0 x 10^^ s; the branching fraction to the two-body decay 
mode is 69.6%. Right panel: The corresponding graph for the antiproton to proton ratio. 
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FIG. 6: Dark matter annihilation diagrams. 

of the Lagrangian that is relevant for computing the rehc density, as well as the dark matter- 
nucleon elastic scattering cross section, is the coupling between x ^"^^ standard model Higgs 

£ D XX^H^H. (4.1) 

In unitary gauge, this can be expanded 

C^^{xhh'' + 2v,^^xhh). (4.2) 

As a consequence of Eq. (14. 2p . x ^i-iid x pairs may annihilate into a variety of standard 
model particles. The leading diagrams are shown in Fig. El The cross section for annihila- 
tions into fermions is given by 



3 



(s — AmlY 

^ ^' (4.3) 



X 

where Nc is the number of fermion colors {Nc = 1 for leptons and Nc = 3 for quarks) and 
is the fermion mass. The cross sections for annihilations into W and Z bosons are given 

14 



by 

^xx^zz = ^7(7^^^2(3 - — + 73t)a/-^' (4-4) 

^xx^w^+H-- - 4^ 3 _ ^2)2 ^2^ + 4m4, ^ y s - 4m2 ' ^^'^^ 

where mw {fnz) is the mass of W {Z) boson. In the case where the dark matter annihilates 
into a pair of standard model Higgs bosons, we can safely ignore the t- and u-channel 
diagrams since the typical momenta are much smaller than at temperatures near freeze 
out. Hence, the cross section is given by 




s-Aml / 6m? 9mf \ 

The evolution of dark matter number density, n^, is governed by the Boltzmann equation 

^ + 3i/(tK = -MK-(nfn, (4.7) 

where H{t) is the Hubble parameter as a function of time and n^'^ is the equilibrium 
number density. The thermally-averaged annihilation cross section, (crv), can be calculated 
by evaluating the integral jsst 

1 /"°° 

<^"> = SmiTKUn,,/T) LJ^'"^ " ^'"xjVSA-.CVS/r) ,s , (4.8) 

where atot is the total annihilation cross section and the Ki are modified Bessel 
order i. We find the freeze out temperature, T/, using the freeze-out condition 

r n^'^{av) 



unctions of 



34| 



1 , (4.9) 



Hitp) H{tp) 

where equilibrium number density as a function of temperature is given by 



rp\ 3/2 

27 



nlQ=('^] e-^^l^ . (4.10) 



The Hubble parameter may be re-expressed as a function of temperature T 

H = l.QQgl^^Tympi. (4.11) 

where is the number of relativistic degrees of freedom and mpi = 1.22 x 10^^ GeV is 
the Planck mass. It is customary to normalize the temperature with the dark matter mass, 
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mpim-^ I —^{av)dx, (4.12) 



X = m^/T. For the points in parameter space discussed below, we found that the freeze 
out happens when Xf ~ 28. The present dark matter density can be calculated using the 
relation 

1 1 [T r gl"' 

where Y is the ratio of number to entropy density and the subscript denotes the present 
time. The ratio of the dark matter relic density to the critical density Pc is given hj Qd = 
2 YoSom^/ Pc where sq is the present entropy density, or equivalently 

VlDh^ ~ 5.6 X 10^ GeV~^ Fq . (4.13) 

Note that the factor of 2 included in the expression for Q^, takes into account the contribution 
from X particles and x antiparticles. 

In the case = 2.5 TeV, we find numerically that the dark matter-Higgs coupling 
A = 0.9 in order that Q^h"^ = 0.1. For = 5 TeV, we find A = 1.8. These order- 
one couplings are perturbative. One should keep in mind that the physics responsible for 
dark matter annihilations is not directly linked to the mechanism that we have proposed to 
account for dark matter decay; other contributions to the total annihilation cross section 
can easily be arranged. For example, if the Higgs sector includes mixing with a gauge singlet 
scalar S such that there is a scalar mass eigenstate near 2m^, then the annihilation through 
the s-channel exchange of this state can lead to a resonantly enhanced annihilation channel, 
as in the model of Ref. 8j. In this case, the correct relic density could be obtained for 
smaller A than the values quoted above. 

Finally, we confirm that the model does not conflict with bounds from searches for dark 
matter-nuclear recoil. In this case, the most relevant contribution comes from the interaction 
between the dark matter and quarks mediated by a t-channel Higgs exchange. The effective 
Lagrangian is given by 

C = -^X^Xqq. (4.14) 



Following Refs 35|, l36||, we can write an effective interaction between the nucleons and dark 
matter, 

^ = -{fpX^XPP + fnX^Xnn) , (4.15) 
where /at = my^/yA/m^, for = p or n. The coefficient An can be evaluated using the 



results of Ref. 



36( 1 : numerically, one finds fp^fn^ ANfnNX/m1 with An ~ 0.35. Given 
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the effective dark matter-nucleon interaction, we find that the spin-independent cross section 
is given by 

47r ml[m^ + rrii^)^ 

For both of the cases discussed earher, (m^ = 2.5 TeV, A = 0.9) and (m^ = 5 TeV, A = 1.8), 
we find asi ~ (9(10~^ cm^. This is two orders of magnitude smaller than the strongest 



bounds, from CDMS 



37|, which range from ~ 2 x 10 cm^ at = 1 TeV to 2 x 10 cm^ 



at m-. 



10 TeV. 



V. CONCLUSIONS 



Models of decaying dark matter require a plausible origin for the higher- dimension oper- 
ators that lead to dark matter decays. The data from cosmic ray experiments like PAMELA 
and Fermi-LAT require that these operators involve lepton fields preferentially. We have 
shown how the desired higher- dimension operators may originate from Planck-suppressed 
couplings between a TeV-scale scalar dark matter particle x vector-like states at a mass 
scale M that is intermediate between the weak and Planck scales. The vector-like sector 
has the structure of a Froggatt-Nielsen model: charged lepton Yukawa couplings arise only 
after these states are integrated out and a discrete gauged Abelian flavor symmetry is bro- 
ken. Couplings between x and the standard model gauge-invariant combination LlHcr are 
then also generated, with coefficients of order (0)^/ (M^ M), where (0) is the scale at which 
the flavor symmetry is broken. Taking M and (</>) near the geometric mean of the reduced 
Planck scale and the weak scale, 0(10^*^) GeV, leads to the desired dark matter lifetime. 
Neutrino masses can be generated via a conventional see-saw mechanism with the mass scale 
of right-handed neutrinos also near M. We pointed out that the symmetry structure of our 
model leads to an overall suppression factor multiplying the charged lepton Yukawa matrix, 
but does not constrain the standard model Yukawa textures otherwise. Hence, our frame- 
work is potentially compatible with a wide range of possible solutions to the more general 
problem of quark and lepton flavor in the standard model. 

We presented the necessary PYTHIA simulations to confirm that our model can account 
for the anomalies observed in the cosmic ray experiments discussed earlier. The leading 
contribution to the primary cosmic ray electron and positron flux in our model comes from 
two-body decays, in which the Higgs field is set equal to its vev in the operator described 



17 



above; the subleading three body decays, x ~^ i~^i~h^, are also possible. We have checked 
that these decay channels do not lead to an observable excess in the spectrum of cosmic 
ray antiprotons, since the cosmic ray antiproton flux is in agreement with astrophysical 
predictions. 

Our model demonstrates that the desired lifetime and decay channels of TeV-scale scalar 
dark matter candidate can be the consequence of renormalizable physics at an intermediate 
lepton flavor scale and gravitational physics at M*. This presents an alternative scenario to 
the one in which dark matter decay is a consequence of physics at a unification scale located 
somewhere between M and M*. 
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